Model of Quantum Criticality in He^ bilayers Adsorbed on Graphite 
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Recent experiments on He^ bilayers adsorbed on graphite have shown striking quantum critical 
properties at the point where the first layer localizes. We model this system with the Anderson lattice 
plus inter-layer Coulomb repulsion in two dimensions. Assuming that quantum critical fluctuations 
come from a vanishing of the effective hybridization, we can reproduce several features of the system, 
including the apparent occurrence of two quantum critical points, the variation of the effective mass 
and coherence temperature with coverage. 

PACS numbers: 71.27.+a, 72.15.Qm, 75.20.Hr, 75. 30. Mb 



He'^ layers, adsorbed on graphite have attracted sub- 
stantial interest since the early eighties, for their remark- 
able properties of surface magnetism [l|, |3, S I3| as well as 
a model system for quantum wetting transitions!^. One 
layer of He'^ atoms adsorbed on graphite pre-plated by 
HD-bilayer or compressed solid has been shown to solid- 
ify at a coverage of 4/7 of the substrate[l, 0, This 
transition has been identified as a Mott transition for 
the He^ fermionsQ. One of the leading questions in this 
field has been to know whether the ground state in the 
Mott phase orders magnetically or is a spin liquid[6|], and, 
in the latter case, which kind of spin liquid -gapless or 
gapfuU- it is0, Theoretical studies have shown that, 
depending on the relative strength of the ring exchange 
parameters[l, [lH, a spin liquid phase, a canted phase, 
or even ferromagnetic ground state can occur. Close to 
the Mott transition, a gapless spin liquid ground state ap- 
pears to be a reasonable choice |12|. Experimental studies 
of the solidification of the first He'^ layer when a second or 
a third layer is adsorbed have been performed long ago 
S There, an enhancement of the static spin sus- 

ceptibility as well as of the specific heat coefficient was 
observed close to the Mott transition. In the experiment 
[H , two layers of He^ are adsorbed on Graphite pre-plated 
by two layers of He^ . The originality of the data T| lies in 
that it's the first system for which, when the second He'^ 
layer arrives at promotion, the first layer is not yet solid- 
ified. As a function of the layer coverage, promotion oc- 
curs at uq ~ 6.3nm^^ while the first layer's solidification 
occurs at ric — 9.9nrn^^ (which corresponds roughly to 
a ratio of 13/19 or 12/19 between the first He^ layer and 
the He'' substrate Q). From specific heat measurements, 
the effective mass is seen to be enhanced like m/m* S, 
with 6 = l — n/uc, and the coherence temperature, below 
which the Fermi liquid behavior is recovered, is shown to 
decrease like Tcoh ~ <5^'^ while approaching the quantum 
critical point (QCP). One striking feature of the data is 
that, strictly speaking, the experiment doesn't reach the 
QCP; at Til ~ 9.2nm^^ the specific heat behavior shows 



hints of a first order transition. Moreover NMR studies 
show that the field-driven magnetization abruptly starts 
to grow at n = 9.2nTO^^. An activation gap extracted 
from the low energy behavior of the specific heat coeffi- 
cient seems to vanish before the quantum critical cover- 
age is reached. In conclusion, the phase diagram of He'^ 
bi-layers pre-plated on two He'' layers, seems to exhibit 
two mysterious phase transitions. One at which the mag- 
netization starts to grow, which corresponds to input of 
a first order transition, and one at which the ratio m/m* 
and the coherence temperature curves extrapolate to zero 
as a function of coverage. 

In this letter we model the system with the Anderson 
lattice in two dimensions with the addition of inter- and 
intra-layer Coulomb repulsion [l4|. Quasi-local f-fermions 
are identified to the first layer He'^ atoms while light c- 
fermions are the second layer He"^ atoms. The bare hy- 
bridization corresponds to the hopping between the two 



layers. In the spirit of the early studies of bulk He^ [15 1 
we solve this model for infinite Coulomb repulsion U be- 
tween the f-fermions, using a slave boson technique equiv- 
alent to Gutzwiller's variational approach. The QCP of 
this system is identified with the Kondo breakdown QCP 
(KB-QCP) of the Anderson latticefH, [13]; that is the 
fixed point for which the effective hybridization vanishes 
and, at the same time, the f-fermions localize. On the dis- 
ordered side of the transition, a spin- liquid phase is neces- 
sary to stabilize the KB-QCP. We can reproduce the vari- 
ations of the inverse effective mass and coherence temper- 
ature with coverage (namely m/m* ~ 5 and Tcoh — 
) in an intermediate temperature regime. Those varia- 
tions are in good agreement with experiments, as seen in 
Fig. 2 where we fit experimental data. While approaching 
the experimental critical coverage, the order parameter 
suddenly drops, reaching the true theoretical QCP before 
the experimental critical coverage is reached. We believe 
this sudden drop of the effective hybridization explains 
the mysterious observation that the system seems to ex- 
hibit two QCPs. Indeed, for us the true QCP is the one 
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at which the hybridization goes to zero, which identi- 
fies experimentally with the point where the magnetiza- 
tion starts to grow at n = 9.2nm~^. The experimental 
QCP corresponds to the extrapolation to T = of the 
inverse effective mass and coherence temperature curves 
obtained in the intermediate energy regime. The drop of 
the order parameter obtained in our model is so abrupt 
that it may trigger as well first order transitions in the 
low energy regime. 

Our starting point is the Anderson lattice model with 
inter- and intra-layer Coulomb repulsion: 



H 



+ 4a iU] - [fLcra + h.C. 

ia 

+ E(^"/.» + ^l"/.*"='* + ^2<i) , (1) 



where (i, j) are the lattice sites created by the Graphite's 
corrugate potential, /^(/) are the creation (annihilation) 
operators for the first layer's fermions, c^(c) arc the 
creation (annihilation) operators for the second layer's 
fermions, ti 
constant, t^. 



: t is the c-fermion's hopping taken as a 
at is the f-fermion's hopping term, V is 
the hybridization corresponding to a hopping term be- 
tween the two layers, i?o < is the f-level potential and 
jjL is the c-fermions chemical potential. Physically, the lo- 
calization of the first layer is driven by the Graphite's cor- 
rugated potential, inducing a triangular lattice in the He** 
bi-layers, inducing itself a second triangular lattice for the 
first He^ layer, commensurate with the substrate's one at 
the "magic" filling number 13/f9[l|. With respect to the 
"13/19" lattice, the f-fermion are thought of being close 
to a Mott transition; the f-band is close to half-filling and 
the hard core Coulomb repulsion leads to strongly corre- 
lated effects. We treat the effect of strong correlations by 
introducing one Coleman's slave boson 18[ which decou- 
ples the /'I' creation operator at site "i" in the following 
way: j}^ fla^i where the f-spinons and the b-holons 
are subject to the constraint J2cr f la fie + ^l^i ~ 1- The 
constraint is taken into account in a Lagrangian formu- 
lation through a Lagrange multiplier A. The properties 
of the second layer c-fermions are very close to the ones 
of the bulk He'^ 1 15j : the Coulomb terms Ui and U2 thus 
mainly re-normalize the hybridization and hopping pa- 
rameters, inducing a dependence in the coverage through 
V = Vo + Vin and t = to + tin, where n is the total cov- 
erage in He^. 

Performing the slave boson decomposition obtains the 



following Lagrangian [l9j : 

^ = E [fl {(^r + Eo + A)% + b.,atb])) 

+ cla ({dr - ^J■)S^J + t) Cj„ 



(2) 




Sj — Uinj/A 



where S — J^ap fl'^fp the spin operator expressed in 
terms of the spinous only, with a the Pauli matrix and dr 
the partial derivative in imaginary time. J = 2{at)'^ /U is 
generated by a second order expansion of our model in U. 
Alternatively, the J-term can be included "ab initio" in 
the model in consideration of the various ring exchange 
parameters generated [1]. We assume that short range 
interactions stabilize a spin liquid with short range ferro- 
magnetic character, in agreement with previous studies. 

In order to fit the experimental data, we evaluate the 
dependance of the bare parameters in coverage. Follow- 
ing 211, we identify Eq as the difference in the average 



potential energy between the first and the second layer. 
Here, the potential energy comes from the joined effects 
of a) Van der Vaals potential between the Graphite sub- 
strate and the layers Vs{z) = {ACl/{27D^)) - 

Ca/z^) where C3 = 2092KA^ is the Van der Waals con- 
stant and D — 192K Q, and b) the Bernardes-Lennard 
Jones [131 potential acting between two He particles 
VLjir) = 4e((cr/r)i2 - (cr/r)6) with e = 10.2K and 
a = 2.56A is the hard core radius. Thus, Eq = En — E12 
with Ell = Vs{zi) + ^^Pi J rdrVLj{r), where i indexes 
the contribution from the various layers(idem for Ei2). 
We find in Kelvin En = -10.495 - 1.07 n{K) while 
E12 — —8.73 — 0.12 n{K) with n the total coverage den- 
sity in nm~^ in good agreement with [2l] |. The param- 
eter Ea then reads Eq = -1.79 - 0.95 n{K). The half- 
bandwidths Dj = 2at = 27r/m/ and Dc = 2t = 2Tr /nic 
are evaluated from Q where the dependence in density is 
extracted from thermodynamic studies of the bulk. We 
find that the dependance in coverage is negligible com- 
pared to the one of Eq; Dc = I.IK, Df — 0.6K thus 
a = 0.55. The relatively high value of a is to be con- 
trasted with the typical values obtained in standard An- 
derson lattice for rare earth compounds. In this system 
the f-band is not particularly flat compared to the c- 
band. The parameter J of the order of a few mK can 
be extracted from the experiments [l[. We take here 
J — AmK. The main difficulty resides in evaluating the 
hybridization V — Vq + Vi n. Noticing that the inter- 
layer spacing is twice smaller than the distance between 
the intra-layer sites, it is reasonable to expect that the 
hybridization is bigger than the c- half bandwidth. In 
this paper we have adjusted the values of Vb and Vi to 
fit the experimental data. We find V = 13.05 — 1.2 n 
, so that at the QCP Vc — 1.89K is larger than Dc, in 
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FIG. 1: Mean-field phase diagram for the Anderson lattice 
model in D = 2 applied to He'^ bi-layers. Following 0] 5 = 
1 — n/ric with ric = 9.9nm~^. The effective hybridization 
Vb drops suddenly at 5 = 0.072, indicating the real QCP. 
The experimental QCP is obtained by extrapolation of Vb 
to T = 0. The Kondo gap A (in red; color online) vanishes 
before the real QCP. 



agreement with the above observation. 

The mean-field equations are obtained by making a 
uniform and static approximation for the hofon cre- 
ation (annihilation) operators b'^{b) in then evalu- 
ating the free energy F and solving for the two equations 
dF/db = and dF/dX = 0. We plotted in Fig.l the 
effective hybridization bV as a function of the coverage 
(the unit for the coverage is identical to the experimen- 
tal ones i5 = {uc — n)/nc with ric = 9.9nm~^ ) and the 
"Kondo gap" A defined as the energy difference between 
the upper band and the Fermi energy. The mean-field 
phase diagrams has two main features. First we observe 
an "elbow" in the order parameter as a function of the 
coverage, corresponding to the set up of the Kondo phase 
(strong hybridization regime). The sharp change of be- 
havior observed corresponds to the emptying of the upper 
band. Note that the opening of the Kondo gap occurs at 
the same point. The model gives an explanation for the 
mysterious observation that the field-induced magnetiza- 
tion (or static spin susceptibility) starts to grow before 
the experimental QCP is reached. In our model, the 
point at which the magnetization starts to grow is the 
physical QCP, which differs from the one obtained ex- 
perimentally, which corresponds to the extrapolation of 
the linear regime to T = 0. The static magnetic suscep- 
tibility is expected to grow quickly in as soon as the first 
layer localizes, since the spin liquid parameter is small 
J ^ AmK. Moreover, from our theory the Kondo gap 
has to vanish before the QCP is reached. This comes 
from the observation that at half filling (corresponding 
to the real QCP), the f-band is half filled, hence con- 
straining the bottom of the first layer to sit below the 
Fermi level. This feature is observed experimentally, if 
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FIG. 2: Coherence temperature (K) , inverse effective mass 
m/m* and effective spinon bandwidth -De// = a'Dc (K) in 
the Anderson lattice model for the He^ bi-layers. a' = b'^a + 
J/t. The dots are experimental data from[l||. The effective 
bandwidth sets the upper temperature of the quantum critical 
regime. The fitting parameters for this model are detailed in 
the text. 



we identify the Kondo gap as the activation gap extracted 
from the thermodynamic measurements of 

The mean-field value for the QCP reads J/t — 
exp [EoD,/V^] 



as m 



17| . We understand that any small 



variation in the value of V with coverage has an expo- 
nential impact on the position of the QCP, justifying our 
option of adjusting the value of V to fit the data. 

We now turn to the fluctuations. A striking observa- 
tion inferred from the experimental data, is the absence 
of quantum critical (QC) regime in temperature. A Curie 
law for the spin susceptibility is observed at very low tem- 
peratures in the localized phase and directly above Tcoh- 
Generically, the upper energy scale of the QC fluctua- 
tions is determined by the first irrelevant operator. In 
our model the formation of the spin liquid and of the 
Fermi liquid are the two mechanisms for quenching the 
entropy. Hence the temperature T* below which the en- 
tropy Rln2 is quenched goes like Max[Def f,Tcoh], with 
Def f the effective bandwidth of the spinous and Tcoh the 
coherence temperature coming from the quantum fluctu- 
ations. The dependence of Dgff and Tcoh with coverage 
are depicted in Fig 2. We see that -De// is of the order 
of J ^ 4mK in the localized phase and follows Tcoh in 
the Fermi liquid phase. The fact that J in this system 
is a remarkably small energy scale, compared to heavy 
fermions, is thus the reason why the QC regime is re- 
duced to much lower temperatures. Quantum fluctua- 
tions are thus most clearly seen through the variation of 
the effective mass and the shape of Tcoh- 

The quantum fluctuations are in the same universality 
class as the ones of the Kondo breakdown model The 
fluctuation spectrum in the intermediate energy regime 
admits the dynamical exponent z = 3, 



q 



a'q 



(3) 
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with Do = 4:kp/{poV'^), ^ is the correlation length as- 
sociated with the fluctuations of b, 7 = mV^ Dq/ (nvp) 
a' — b'^a + J/t and po — mc/(27r) is the c-fermions den- 
sity of states. The boson mass is given by nib = D^^ 
evaluated at T = 0. We evaluate it by differentiat- 
ing twice the mean-field equations with respect to the 
bosonic field b and evaluating the result at the mean- 
field saddle point. For the effective mass, we use the 
Luttinger-Ward expression of the free energy typical of 
Heisenberg-type theories, analogous to the one derived 
'^F = Ffg + T/2 S„ / d2q/(2^)2 log [D-\q, 



where Fpo is the free-energy of the system at the mean- 
field. Since the system at the mean-field consists of two 
hybridized bands, one obtains after diagonalization 



F = -- 



ttT' 
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(4) 



where pi{p2) are the density of states of the upper (lower) 
bands. The effective mass thus reads m* = 27r(pi -I- P2) + 
7^/(4q!'). The coherence temperature is computed by 
evaluating the corrections to scaling to the boson propa- 
gator ([31); namely the one loop diagrams responsible for 
the temperature dependence of mi,{T) — D^^ ^~^(T). 
Following [2^ we find (remember z — 3 here ) 

mf,(T) = jUbiT = 0) + C TLogT , (5) 

where C had to be adjusted to C = 7.5 10^'^ to fit the 
data, while the analytic evaluation gives C = DJ/{6V^). 
The coherence temperature Tcoh obtains when the cross- 
over condition mb{T) = is satisfied. The results for 
m/m* and TcoH are presented in Fig. 2 and directly com- 
pared to experiments. Since we work within a slave boson 
saddle point approximation, the need of one fitting pa- 
rameter for the amplitude of the coherence temperature 
is to be expected. The exponents can be understood in a 
simple way. For z = Z theories in the Fermi liquid phase. 
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the effective mass goes like the correlation length 
m/m* ~ ^ 

we see that £, " ^ \Jmi, 
ature goes like b^. In the regime where b varies linearly 
with the coverage n we thus get 



^. From the dispersion of the boson mode 
1 . /™. - \) Now the coherence temper- 



m/m* 



est - 



Tr. 



{est 



(6) 



In conclusion we have performed a study of quantum 
criticality in He^ bi-layers by mapping the experiments 
on an extended version of the Anderson lattice in two 
dimensions. We examined the possibility for the Kondo 
breakdown QCP to be responsible for the quantum fluc- 
tuations observed. The system of He^ bi-layers enables 
to directly test the theory from the bare parameters. Our 
model is successful in that 

• it explains the occurrence of two QCPs seemingly 
observed experimentally. Our interpretation is that 
one of the QCPs observed experimentally corre- 
sponds to the extrapolation to T = of an inter- 
mediate energy regime; 



it predicts that the activation gap extracted from 
thermodynamic measurements vanishes before the 
QCP is reached; 



• it gives exponents for the variation of the effective 
mass m/m* and the coherence temperature Tcoh h^ 
good agreement with experiment, so that a direct 
fitting of the data is possible. 



We have used three fitting parameters to account for 
the prefactors of the effective mass and coherence tem- 
perature as well as to precisely determine the position of 
the QCP from the bare parameters. Our study is the first 
case where an itinerant QCP showing non Fermi liquid 
behavior is used to fit experimental data from the bare 
parameters. 
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